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Q \ Abstract 

It is shown that a scalar field, minimally coupled to gravity may have collapsing modes 
even when the energy condition is violated, that is, for [p + 2>p) < 0. This result may 
be useful in the investigation of the possible clustering of dark energy. All the examples 
dealt with have apparent horizons which form before the formation of the singularity. The 
singularities formed are shell focusing in nature. The density of the scalar field distribution is 
seen to diverge at singularity. The Ricci scalar also diverges at the singularity. The interior 
spherically symmetric metric is matched with exterior Vaidya metric at the hypersurface 
^ ■ and the appropriate junction conditions are obtained. 
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c<i '. 1 Introduction 



The ultimate fate of a collapsing matter distribution remains a subject of fervent interest for 
^ ■ many years. It is generally believed that a collapsing matter distribution, if it satisfies the energy 
condition (p + 3p) > 0, will crush to a singularity. There is however no general indication whether 
an event horizon would form to shield this singular state of matter from a distant observer. In 
fact, there are ample theoretical examples where the singularities could be naked to an external 
observer either for a finite period of time or even for ever. For a comprehensive review we refer 
to the monograph by Joshi[l]. A lot of work therein or that carried out later normally considers 
the collapse of a perfect fluid. As the final stages of collapse might involve dissipative processes 
as well, some later investigations also consider viscous effects (see [2] and references therein). 

Investigations on collapsing gravitational systems with a scalar field revealed that the non 
linearity of Einstien equations could lead to critical phenomena close to the threshold of black 
hole formations [3]. Naturally this led to an arena of great interest (for recent reviews we refer to 
[4]). Furthermore, the scalar field collapse could also lead to the formation of naked singularities. 
This was demonstrated by analytical investigations by Christodoulou[5], Goswami and Joshi[6], 
Giambo[7] and also by numerical work of Choptuik[3], Brady[8], Gundlach[9], and others. 

In a series of papers Christodoulou[10, 11, 12] pioneered analytical studies in gravitational 
collapse of a spherically symmetric configuration of massless scalar field models. Detailed inves- 
tigations on such models have been carried out later [3, 8, 13]. Collapsing spherically symmetric 



^E-mail: koyel_g_m@yahoo.co.in 
^E-mail: narayan@iiserkol.ac.iii 



2 



massive scalar field models have also appeared in the literature[7, 14, 15, 16]. Gonclaves[14] in- 
vestigated the black hole formation with a massive scalar field in an Einstein-de Sitter universe. 
Conclaves and Moss[15] also showed that for a range of parameters, the collapse of a spherically 
symmetric scalar field can be formally treated as a collapsing dust ball. This is similar to the 
case of 7 = I in the present work (section 3.3). Giambo[7] showed that a massive scalar field, 
if collapses completely, can give rise to a naked singularity. Goswami and Joshi[6] constructed a 
class of solutions with a massive scalar field which could lead to a naked singularity. 

The very recent interest in scalar field collapse stems from a cosmological requirement. It 
has already become a folklore that the present Universe is undergoing an accelerated phase of 
expansion. Scalar fields endowed with a potential are amongst the most favoured candidates as the 
"dark energy"- the agent driving the cosmic acceleration [17]. Although the general speculation 
is that the dark energy does not cluster at any scale below the Hubble scale, there is not much 
evidence for this. The collapse dynamics of a minimally coupled scalar field with a negative 
effective pressure can indeed help us to have a better understanding of the scenario. 

The objective of the present investigation is to explore the collapsing modes of a simple 
spherically symmetric minimally coupled scalar field model. The scalar field is assumed to 
depend only on time, so this is an analogue of Oppenheimer - Snyder model of the collapse of a 
spatially homogeneous fluid [18]. 

Section 2 describes the scalar fleld distribution and obtains the general flrst integrals in sit- 
uations where the effective pressure and energy density of the scalar field are related as 
Pff, — wp^, w being a constant. Section 3 discusses special cases of the collapsing models. In 
section 4 the collapsing spacetime is matched at the boundary with an external metric and in 
section 5 some concluding remarks are made. 

2 Spherically symmetric time dependent distribution of a 
scalar field 

The relevant action for a scalar field minimally coupled to gravity is given by 

A^l V^d^'xlR + ^cf>^y,^ - VicP)]. (1) 

We take the spherically symmetric metric as 

ds^ = dt^ - A\r, t)dr^ - B\r, t)dn^, (2) 

where A and B are functions of the radial coordinate 'r' and time 't'. A minimally coupled scalar 
field endowed with a potential V — V{(p) is considered. For the sake of simphcity, the scalar 
field is taken to be spatially homogeneous and evolving with time alone. The non-trivial field 
equations become 

2B" 2A'B' B'^ Ab B^ 1 1 

+ ^n^-^ + 2-r^ + ^ + ^ = o'/'' + n0), (3) 



A'^B A^B A'^B^ AB B'^ B^ 2 
B'^ B B^ 1 1 



2o-^-^ = o'^ -W). (4) 



^2^2 B B^ B^ 2 
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B" A'B' A B AB l-^ 

(5) 



A^B A^B A B AB 2 
B' A B' ^ 

where (j) = (f)(t) is the scalar field, V{(f)) is the scalar potential. The equations are written in units 
where SttG = 1. An overhead dot and prime indicates difi^erentiation w.r.t 't' and 'r' respectively. 
As already mentioned, the contribution from matter (the scalar field in this case) is spatially 
homogeneous like the homogeneous fiuid collapse discussed by Oppenheimer and Snyder [18]. 
Equation (6) readily integrates to yield 

B' = X{r)A, (7) 

where A(r) is an arbitrary function of 'r'. In what follows A(r) is taken to be constant for the sake 
of simplicity. This indeed is a special case, but these would lead to some tractable solutions so 
that the possibihty of collapse can be investigated. However, if A{r, t) is separable as a product 
of functions of V and , the assumption (7) is hardly restrictive as A(r) in that case can be 
absorbed by the scaling of the radial co-ordinate. Anyway, the physical motivation is that in this 
case the collapsing model will eventually become an FRW one. If A and B are now written as 
A"^ = e" and -B^ = e^, the field equations take a much simpler form as 

f + + (8) 

ii + li^' = -l<t>' + Vi-P). (9) 
2+T + f + T + T = V '^o' 

where equation (7) has been used, which now reads as 

e" = ^^e^. (11) 
4 

The wave equation for the scalar field is given by 

r ;^ + V^ = Q. (12) 

In what follows, the potential V{(f)) will be assumed to be proportional to 0^ , 

^(0) = 70', (13) 

where 7 is a constant. Indeed this choice does not have any physical requirement, but along with 
a major simplification for the equation system the ansatz provides the fourth equation which we 
require to solve the system of equation. We have only three independent equations (8)- (10) to 
solve for four unknowns a, P, (f> and V. As goo — 1 and = 0(^), 0^ is actually which is 

indeed a scalar and equation (13) is thus consistent. One can see that 7 > | would lead to the 
case where the effective pressure is negative. 

However, it deserves mention at this stage that the field in this case is not a k-essence field 
where the scalar field Lagrangian is given by L{X), where X = (V0)^. Actually the model does 
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have a potential and the present work deals with those which are proportional to the kinetic part. 
We shall see later that an exponential potential is in fact consistent with the ansatz(13) in the 
present model. 

The equation of state parameter, w — ^ — i^2^|^|^| will be a constant given hy w = j^f^- 

For 7 > 1, the equation of state parameter u is more negative than — | and the scalar field acts as 
a dark energy. For any finite positive value of 7, w remains greater than —1 and hence the scalar 
field is in fact a quintessence field and cannot give a phantom field with w < —1. If 7 = 0, that 
is, for the case without potential, the value of w becomes 1 and the field is then formally similar 
to a stiff fluid. For 7 = |, one has an effectively zero pressure, and the scalar field is formally 
equivalent to a pressureless dust. 

With the help of equation (13), the wave equation (12) yields a first integral as, 

- Jr^sir), (14) 

where g{r) is an arbitrary function of r to be determined from initial conditions. As (p is assumed 
to be spatially homogeneous, this equation indicates, 

e(f+/3)^^(^)/(i), (15) 

where f{t), a function of time, and has to be determined from the field equations. This gives the 
form of equation (14) as, 

<i>^f{t)-^. (16) 

Using (15) in a combination of field equations (8)-(10) one obtains 

f-umf^, (17) 

which readily yields a first integral as 

/2 = ^(27 + l)/^+Ci, (18) 

Ci being an arbitrary constant. In order to satisfy the field equations, one can verify that Ci = 0. 
Choosing g{r) as the solution of e" and are given as, 

e« = (/(^))2/^ (19) 

and 

= mr^'r'. (20) 
Equation (18) can be written as (with Ci — 0), 



/ = -y^(27 + l)/^- (21) 

The negative root is chosen in order to consider a collapsing scenario. Equation (21) readily 
integrates to yield 
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where to is another constant of integration. 

The solution of the scalar field as some function of t is given by (using equations (14) and 
(15)), ^ 

= -\^^ln(^ - ,r^t) . (23) 
^ V 3 V27 + 1 V 2(27 + 1) ; 

As the complete set of solutions is now in hand, it is not difficult to check which form of 
potential is consistent with the ansatz(13). It is easy to see from equation (23) that 
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V{4,)=^exp{ J/,„ |. (25) 



In conjunction with equation (13), one can now write, 

3 Collapsing models 

3.1 7<l/2 

For values of 7 < 1/2, the matter content satisfies the energy condition (p + 3p > 0). Choosing 
7 = 1/4, one has the solution of f{t) as 

/2 

fit) = (3^0 - t) (26) 

Now et+^, which according to equation (15) is given by g{r)f{t), determines the proper volume. 
The above equation clearly shows that f{t) indeed has a non zero value for all t less that ^ but 
becomes zero when time t = |to- So for this value of t (designated as tg), one has a singularity 
of a zero volume. Choosing a value of tg — 3, one gets tg — 2. The plot of f{t) vs t in Figure 1 
clearly reveals the collapsing scenario. 

The condition for the apparent horizon is 

g'^'^RaR^ = 0, (27) 

where R is the proper radius of the two sphere. In the present case R — B — e^^'^ and equation 
(27) reads like 

^oos2 + ^iiS'2^0, (28) 

that is 

^2 _ ^ ^29) 

which in view of (11) yields 

= 4. (30) 

The time of formation of the apparent horizon , tah, has been calculated for a fixed value of r = 1 
in all the cases discussed here. In this case (7 = 1/4 and to is chosen as 3), tah = 1-967. Hence 
the singularity formed is shielded by the horizon. 



6 



3 

2.5 
A 2 
1 1.5 
1 1 

0.5 



0.5 1 1.5 2 

>t 

Figure 1: Plot of f{t) vs. i for 7 = 1/4 putting to = 3. 

3.2 7 > I 

Choosing a value of 7 > | makes p + 3p < 0, that is the energy condition is violated and the 
scalar field is apt to act as a dark energy. With a choice of 7 = 3/2 (as an example), the solution 
of f{t) turns out to be, 




It is evident from the above expression that f{t) becomes zero when t = Hence one indeed 
has a singularity of a zero proper volume (/ = 0) at a finite future. Choosing to = 4 one gets 
ts = 1.632 and tah = 0.636. Figure 2 shows the plot of f{t) against time for choice of to = 4. It is 
observed that initially there is a non zero proper volume and then the singularity of zero volume 
is attained at a finite time. A change in initial condition changes the time at which singularity 
is attained keeping the nature of the plots unchanged. The time at which apparent horizon is 
formed also changes with the change of initial condition. For the choice of to=4, tah < tg- Hence 
the singularity formed is covered by the horizon from an external observer. 
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Figure 2: Plot of f{t) vs. t for 7 = 3/2 putting to = 4. 

For some other choice of 7 > 1/2, for example 7 = 5/2 and choosing to — 6, f{t) vs t is plotted 
as shown in Figure 3. The plot shows that in this case also f{t) becomes zero at some finite time 
(tg = 2). The value of tah is 0.592. In all the cases tah < ts, thus the singularities formed are 
hidden behind the horizon. 
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Figure 3: Plot of f{t) vs. t for 7 = 5/2 putting to = 6. 



The energy density and the effective pressure of the scalar field are given hj = (| + 

l)f(t) ^'^''^'^ s-iid = {2 ~ 7)/(^) respectively. From the expressions for f(t) it is clear that 
for any positive value of 7 whether it is greater than or less than 1/2, the density and pressure of 
the scalar field would diverge at the time at which the proper volume becomes zero. The Ricci 
scalar R also diverges at the singularity as is evident from its expression given below, 



f 9P ■ 



(32) 



3.3 7 = 1/2 

For 7 = 1/2, the model formally resembles a perfect fluid distribution with zero pressure. From 
equation (22), one has the solution of f{t) as 



{to-Vsty 



(33) 



At time t = ^ the proper volume becomes zero. Choosing to — 2, one gets ts = 1.155. Fig 4 
shows the plot oi f{t) vs t for 7 = 1/2. The time of apparent horizon formation in this case is 




1.2 



Figure 4: Plot of f{t) vs. i for 7 = 1/2 putting to = 2. 
tah = 1.071 which is once again less than the value of tg. 
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Singularities formed in collapse can be shell focussing or shell crossing in nature[19, 20]. For 
a spherically symmetric collapse the shell focusing singularity occurs at gee = and the shell 
crossing singularity occurs when = 0, gee and r > 0. It is evident from (2) and (20) that 

996 = f'^fiff^^- At singularity since fit) = 0, therefore gee = at singularity, clearly indicating 
that the singularities formed in all the cases are shell focusing in nature. 



4 Matching with the exterior metric 

The interior metric now has to be matched with an exterior metric at the boundary E. Following 
ref [6], we match the spherical ball of collapsing scalar field to a generahzed Vaidya metric . The 
metric just inside E is, 

dsJ = df - e'^dr^ - e^dO^ (34) 
and the metric in the exterior of E is given by 

ds^^ = (1 _ ^MhA)dv^ + 2dvdr, - r,^dn\ (35) 

Matching the first fundamental form on the hypersurface we get 

' (36) 



(its J I 2M(r„,ii) ^ 2dr^ 

and 

r. = rfitf (37) 
Matching the second fundamental form yields, 

_ 2M{rv,v) I 2dr\, \ 
^ I r ^ do (3g) 

Vv dv ' 



Using equations 36, 37 and 38 one can write 

dv 1-^ 



dt 1 _ 2M/-V3 ■ 
r 

From equation 38 one obtains 



(39) 



^ ^ 9 ^ V ^27273 ^ ^171073 9/2 ^^^^ 

j,2 jljZ 

Matching the second fundamental form we can also write the derivative of M(i>, r„) as 

M 2r^ 

^(^-^U = ;7T7^-9^- (41) 

Equations 37, 39, 40 and 41 completely specify the matching conditions at the boundary of the 
collapsing scalar field. 
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5 Discussions 

The present work investigates a simple spherically symmetric spacetime with a scalar field. The 
scalar field is self- interacting, i.e., it has a potential. The field is chosen as spatially homogeneous, 
so it represents the scalar field analogue of Oppenheimer-Snyder collapse for a perfect fluid [18]. 
At least for the special case where the potential could be written as proportional to the kinetic 
part of 0, it could be shown that the model possesses a collapsing solution. 

Most of the models on massive scalar field collapse, available in the literature [7, 14, 15], 
demonstrate the formation of either a naked singularity or a black hole are subject to the initial 
conditions. They generally work it out for energy condition (p + 3p) > 0. In our work we 
demonstrate the collapse of a scalar field endowed with a potential for both the energy conditions 
(p + 3p) > and (p + 3p) < 0. The value of the equation of state parameter -u;, determined by the 
value of 7, dictates whether the energy condition is satisfied or not. When 7 is less that —1, the 
energy condition is violated. For (p+3p) > 0, like many other models[6, 7, 16], we get singularities 
at a finite time and the possibility of the formation of an apparent horizon ahead of the collapse 
to the singularity. But the intriguing feature of the present work is that we obtain singularities 
at a finite time even for the condition (p + 3p) < wherein the effective pressure is negative. 
So this is an investigation of the collapse scenario of dark energy itself rather than collapse of a 
fiuid against a dark energy background [22]. This is an indication towards the possibility of the 
clustering of dark energy, where (p + 3p) < 0. 

All the examples discussed in the present work allow the formation of an apparent horizon 
well before the distribution crushes to the singularity of zero volume. So the singularities are not 
naked. The interior metric has been matched with the exterior Vaidya metric and the appropriate 
matching conditions are obtained. 

In all the cases (for different values of 7 chosen) it is observed that the energy density as well 
as the effective pressure for the scalar field diverge at the singularity. It has been checked that 
the scalar curvature R = also diverges as the proper volume becomes zero. The scalar field (f) 
also is singular at tg. Equations (26) and (31) have a deceiving feature, as if although the volume 
is zero at t = tg, it becomes quite well behaved beyond tg. A look at equation (23) reveals that 
actually beyond tg the scalar field itself becomes undefined as the argument of logarithm becomes 
negative. So one indeed has a singularity at a finite future, beyond which the model cannot be 
discussed any further. 

The present work deals with a particular potential and docs not represent a general kind of 
potential. Although we started with an assumption that the potential is proportional to the 
kinetic part, actually the model is that of an exponential potential (25). 
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